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Abstract 


Flooding is a natural phenomenon and the risk of its occurrance in urban areas is a global issue. 
Rapid urbanization and climate changes have increased the risk of urban flooding, leading to 
massive infrastructure and human losses. The aim of this study was geomorphological zoning of 
Shahinshahr River flood risk. Therefore, the HEC-RAS numerical model was used to simulate 
the flood. The geometric data was processed in GIS by HEC-GEORAS extension. The return 
periods of 25, 50, and 100 years of catchment area as well as its physiographic characteristics 
including area, length of main stream, CN curve number, concentration time, latency in the 
watershed were entered into the HEC-HMS software. Accordingly, the maximum flood 
discharge with different return periods were calculated whereby the river area was devided into 
three reaches. Finally, manning's roughness coefficient was calculated for each reach that the 
Manning roughness coefficient was calculated using the Coon method. The results show that the 
flood zone in the 25-year return period includes 0.948 km?, 1.13 km? in the 50-year return 
period, and 1.34 km? of the lands along the Shahinshahr River in the 100-year return period. In 
the last reach, due to the reduction of the slope and the flow velocity, the flood zone has become 
wider than other periods. 
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مقاله پژوهشی 
تحلیل دینامیکی مخاطره سیلاب با استفاده از مدل هیدرولیکی HEC-RAS‏ 


(مطالعه مو ردی: رودخانه شاهین شهر استان اصنهان) 
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حکیده 


وقوع سیل یک پدیده طبیعی است و خطر وقوع آن در مناطق شهری یک مسئله جهانی است. رشد 
شهرنشینی و تغیبرات آب و هوایی مخاطره سیلاب‌های شهری را افزايش داده و باعث ایجاد خسارت به 
زیرساخت‌ها و صدمات gl‏ شده است. هدف این پژوهش تحلیل دینامیکی مخاطره سیلاب رودخانه 
شاهین‌شسهر (استان اصفهان) است؛ بنابراین جهت شبیه‌سازی مخاطره سیلاب از مدل یک 
بعدی HEC, RAS.‏ استفاده گردید و از طریق الحاق HEC-GEORAS‏ در محیط GIS‏ پردازش 
داده‌های ژئومتری انجام شد. جهت محاسبه دبی برای دوره بازگشت‌های ۲۵ ۵۰ و ۱۰۰ alle‏ حوضه آبریز 
منطقه موردمطالعه به داخل محیط نرم‌افزار HEC- HMS‏ وارد شد و مشخصات فیزیوگرافی حوضه آبریز 
از قبیل: مساحت. طول آبراهه اصلی. شماره منحنی خاک (CN)‏ زمان تمرکز و زمان تأخیر وارد گردید. بر 
این اساس. نتایج خروجی حداکثر دبی سیلاب با دوره‌های بازگشت‌های متفاوت محاسبه شد نتایج حاصل 
از مدل نشان می‌دهد که با توجه به شرایط رودخانه. منطقه به سه بازه مختلف تقسیم‌بندی شد که ضریب 
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حاصل از مطالعات OL‏ می‌دهد که پهنه سیلاب در دور بازگشت Yo‏ ساله ۰/۹۶۸ کیلومتر مربع» در دوره 
بازگشت ٩۰‏ ساله ۱/۱۳ و در دوره بازگشت ۱۰۰ ساله ۱/۳۶ کیلومتر مربع از اراضی حاشیه رودخانه 
شاهین‌شهر را در بر گرفته‌اند؛ و در بازه آخر به دلیل کاهش میزان شیب و کاهش سرعت جریان پهنه 
سیلاب از سایر بازه‌ها گسترده‌تر شده است. 


کلیدواژه‌ها: شبیه‌سازی سیلاب. رودخانه شاهین‌شهر -HEC_RAS‏ 


۱- مقدمه 


بیشتر جوامع بشری در کنار رودخانه‌های بزرگ جهان مستقر و شکوفا شده‌اند؛ زیرا رودخانه‌ها و سیلاب دشت‌ها 
برای زندگی بشر و به‌ویژه برای کشاورزی و توسعه شهری مهم هستند ( یوسفی " و همکاران. ۲۰۲۱). اغلب 
سیستم‌های رودخانه‌ای تحت تأثیر متفیرهای خارجی ازجمله تغییرات اقلیمی و زمین‌شناسی دچار آشفتگی می‌شوند. 
الگوها و میزان تعدیل‌های رودخانه‌ای و یا تغییرپذیری آن‌هاء ماهیت رویدادهای آسفتگی را نان می‌دهند. این 
آشفتگی مانند یک سیل فصلی می‌تواند پدیدهایی دوره‌ای با وقوع کې شدت زیاد و کوتاه‌مدت باشد و اثر موضعی بر 
رودخانه بر جای بگذارد (حسین‌زاده و همکاران. ۱۳۹۶). سیل یکی از مخرب‌ترین مخاطرات طبیعی؛ بخصوص در 
مناطق حشک و نیمه‌حشک L)‏ سیل‌های ناگهانی) که بارندگی زیادی را تجربه می‌کنند می‌باشد. این مخاطره ناشی از 
پیوند چندین عامل طبیعی و انسانی است. سیلاب هنگامی رخ می‌دهد که دبی رودخانه به حدی افزایش می‌یابد که 
آب از سواحل رودخانه عبور کرده و مناطق مجاور را غرقاب می‌کند. درنتیجه جوامعی را که در مجاورت رودخانه‌ها 
ساکن هستند. تحت تأثیر قرار می‌دهد علاوه بر آن تغییرات آب و هوایی احتمال وقوع سیل و بزرگی آن‌ها را افزایش 
می‌دهد ( ترمه " و همکاران» ۲۰۱۸؛ ازین " و همکاران. ۲۰۲۰). سیل سالانه به‌طور متوسط ۷۵ میلیون نفر را در سراسر 
جهان تحت تأثیر قرار می‌دهد (حدود ۲۰۰۰۰ تلفات در هر سال). ضرر اقتصادی حاصل از آن ۹۵ میلیارد دلار بوده 
است در سال‌های اخیر» فراوانی سیل بیش از 2٤٩‏ افزايش یافته است؛ زیرا تغیبرات آب و هوایی به‌طور قابل‌توجهی 
بر شدت الگو و بزرگای سیلاب‌ها تأثیر گذاشته است. به‌ویژه کشورهای آسیایی با برخورد نسبتاً زیاد در معرض سیل 
و زلزله قرار دارند در آسیا بیش از ۸٩۰‏ تلفات انسانی ناشی از بلایای طبیعی به علت جاری شدن سیل است (چانگ 
و چن Y Y M‏ حسروی" و همکاران؛ ۲۰۱۸ ۲۰۲۱). ایران اخیراً سیلاب‌های فاجعه باری را تجربه کرده است که 
علاوه بر تأثیر تغییرات آب و هوایی. تأثیر فعالیت‌های انسانی ناشی از مدیریت ضعیف حوضه آبریز ازجمله 


جنگل‌زدایی» چرای بیش ازحد دام توسعه سکونتگاه‌های روستایی و شهری باعث کاهش فضای هیدرولیک 
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رودخانه‌ها شده است» سیل‌های مارس و آوریل ۲۰۱۹ در ایران که Yo‏ استان از ۳۱ استان کشور را تحت تأثیر قرار 
داده است. درنتیجه بیش از ۷۷ نفر OWE‏ انسانی و خساراتی به ارزش ۲/۲ میلیارد دلار به سکونتگاه‌های شهر و 
روستایی و زمین‌های کشاورزی وارد شده است -(https://en.wikipedia.org/wiki/2019_Iran_floods)‏ در Ol»!‏ 
حدود ۱۲ درصد اراضی مستعد ابتلا به مخاطره سیلاب هستند و ۳ درصد دارای احتمال خطر سیلاب با ریسک بسیار 
زیاد می‌باشد ( خسروی و همکاران» ۲۰۲۱). در سال‌های etl‏ ارزیابی‌ها و برآورد سیلاب به دلیل افزایش شهرنشینی 
و تمرکز سکونتگاه‌های مسکونی در اطراف رودخانه‌ها و دشت‌های سیلاب به‌طور قابل‌توجهی افزایش یافته است ( 
آوند ' و همکاران» ۲۰۲۱). ارزیابی و مدیریت ریسک سیل یک اقدام مهم برای شناسایی احتمال حطر مناطق مستعد 
سیل و کاهش حوادث احتمالی سیل است. اگرچه ساخت سازه‌های هیدرولیکی نقش مهمی در JES‏ سیل دارد اما 
برای کاهش حسارت ایجاد شده کافی به نظر نمی‌رسد؛ بنابراین» ایجاد پایگاه داده‌ای از نقشه‌های خطر سیل ضروری 
است. سرویس حفاظت از منابع طبیعی آمریکا (USDA)‏ که قبلاً خدمات حفاظت از خاک (SCS)‏ نامیده می‌شد 
روش شماره منحنی خاک (CN)‏ را جهت محاسبه میزان رواناب حوضههای cy pl‏ توسعه داد شماره منحنی 
خاک (CN)‏ از تجزیه‌وتحلیل تجربی رواناب حوضه‌های آبریز و با فاکتورهایی از قبیل شیب کاربری و گروه‌های 
هیدرولوژیکی خاک منطقه محاسبه می‌شود (خدمات حفاظت خاک امریکا ۲ ۱۹۷۲ oe‏ ۲۰۱۸). تحقیقات مختلف با 
استفاده از مدل‌های هیدرودینامیکی سیلاب را در دشت‌های سیلابی شبیه‌سازی می‌کند این مدل‌های عددی» معادئات 
جریان را به‌عصورت یک بعدی و دو بعدی و چند بعدی حل می‌کنند. در میان مدل‌های عددی (HEC-RAS)‏ 
گسترده‌ترین مدل یک بعدی شبیه‌سازی احتمال خطر سیلاب می‌باشد که به‌وسیله برونر (۲۰۰۱) گسترش پیدا کرد و 
جریان‌های هیدرولیک را Col = Cy pee‏ ماندگار و غیرماندگان محاسبات انتقال رسوب و مدل‌های درجه 
حرارت را محاسبه می‌کند مطالعات متعددی با استفاده از نرم‌افزار HEC-RAS‏ انجام شده است که شامل: محاسبه 
سیلاب (آبدررزاک" و Ge ۲۰۰۹ Ohl Saa‏ ۲۰۱۸؛ زلناکوا* و iab Sa‏ ۲۰۱۹؛ گراوند" و همکاران ۲۰۲۰ 
مدل‌سازی انتقال رسوب و حمل بار معلق (پیترو و همکاران. ۲۰۱۵ و بررسی تغییرات مورفولوژی سنگ‌بستر و 
کانال (کلارک" و aL Kas‏ ۲۰۰۸) بوده است. اخیراً GIS‏ به دلیل توانایی OF‏ در مدیریت مقدار زیادی از اطلاعات به 
یک ابزار اساسی در هیدرولوژی تبدیل شده است؛ زیرا این نرم‌افزار توانایی زیادی در مرتب‌سازی مقدار زیادی از 
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از منابع مختلف را دارا می‌باشد و gal‏ 618 با مدل یک بعدی HEC-RAS‏ پهنه گسترش سیلاب را در رودخانه‌های 
مختلف محاسبه می‌کند ( عز و همکاران ۲۰۱۸)؛ و علاوه بر آن تحقیقات زیادی توسط محققان داخلی در زمینه 
مدل‌سازی سیلاب انجام شده است که عبارت است از: دامادی و همکاران(۱۳۹۹) به پهنه‌بندی سیلاب با استفاده از 
مدل هیدرولیکی (HEC-RAS)‏ در رودخانه سرباز استان سیستان و بلوچستان به مطالعه پرداعته است و علاوه بر آن 
روستایی و همکاران (۱۳۹۹) به شبیه‌سازی شبه دو بعدی جریان سیلاب رودخانه لیقوان با تأکید بر دشت سیلابی را 
موردمطالعه قرار داده‌اند. همچنین مختاری و همکاران(۱۶۰۰) به تحلیل دینامیکی مخاطره سیلاب با استفاده از مدل 
هیدرودینامیکی HEC-RAS‏ و تکنیک GIS‏ مخروط افکنه لیلان. شمال غرب ایران را موردمطالعه قرار داده اسست و 
همچنین ریاحی مدوار و فکوری (VE)‏ به تحلیل عدم قطعیت نتایج مدل HEC-RAS‏ در شبیه‌سازی فراسنجه‌ای 
هیدرولیکی جریان رودخانه کارون با رویکرد مونت‌کارلو به مطالعه پرداخته‌اند که در تمامی این مطالعات پهنه سیلاب 
در دوره بازگشت‌های مختلف و میزان آب‌گرفتگی در هر کاربری را محاسبه کرده‌اند. اکوسیستم‌های خشک و 
نیمه‌عشک معمولاً دارای سیلاب شدیدتری نسبت به سایر اقلیم‌ها هستند منطقه موردمطالعه حوضه آبریز شاهین‌شهر 
بوده است شرایط اقلیمی خشک و سرد نوع سازندهای زمین‌شناسی» خاک فرسایش پذیر (متشکل از سیلت و ماسه) 
و شیب توپوگرافی از ویژگی‌های بارز حوضه ابریز شاهین‌شهر بوده است که سبب بالا بودن پتانسیل ALS‏ رواناب و 


درنتیجه فرسایش LAS‏ خاک و تولید رسوب در حوضه مذکور شده است. 


۲- مواد و روش 
= منطقهٌپژوهش 


منطقه موردمطالعه حوضه آبریز شاهین‌شهر دارای مساحت OWE‏ کیلومتر مربع در مرکز ایران و در محدوده 
سیاسی استان اصفهان (در غرب شهر اصفهان) قرار گرفته است. حوضه آبریز شاهین‌شهر شامل حوضه‌های آبریز 
میمه» مورچه‌خورت. اصفهان- برخوار و حوضه آبریز علویجه -دهق که زیر حوضه آبریز گاوخونی می‌باشد این 
حوضه در مختصات ٥۰ TV.‏ تا Y‏ ۵۱۳ طول شرقی و ۳۲۳,۵۲ تا ۳۳۶۱ عرض شمالی واقع شده است. 
حداکثر ارتفاع حوضه از سطح دریا ۳۶۵۸ متر و حداقل ارتفاع نیز MAY.‏ متر می‌باشد. خروجی آب از این حوضه به 
حوضه ابریز کوهپایه - سگزی و درنهایت به تالاب گاوخونی منتهی می‌گردد. حوضه آبریز شاهین‌شسهر از نظر 
ژئومورفولوژیکی در واحد سنندح-سیرجان قرار دارد و رسوبات منطقه متعلق به دوران دوم و کواترنر می‌باشد 
رسوبات دوران دوم شامل دوره‌های تریاس و کرتاسه است که اکثر آن در ارتفاعات منطقه وجود دارد. مسیل 
شاهین‌شهر (بازه موردهطالعه) با طول A‏ کیلومت آبراهه اصلی حوضه آبریز مسیل شاهین‌شهر می‌باشد در پیرامون 
مسیل شاهین‌شهر اغلب شامل نهشته‌های گراولی کواترنری؛ ته نشست‌های رودخانه‌ای (سیلت و ماسه) تشکیل شده 
است. متوسط بارش حوضه در طی دوره آماری ۲۰ ساله ۱۱2/۳ میلیمتر و میانگین (cles‏ ماهانه 3 ole‏ ۱۵/۸ درجه 


سال يازدهم تحلیل دینامیکی مخاطره سیلاب با استفاده از .... ^Y‏ 


سانتی کر اد در کل محدوده موردمطالعه می‌باشد. در حوضه از اوایل آبان تا اواسط اردیبهشت مرطوب بوده و سایر 


aul 02 9 نفر‎ YE£VW از نظر جمعیتی محدوده شاهین‌شهر و ميمه دارای جمعیتی بالغ‌بر‎ Yao 


Morse Lus uci Sr'2s25e 
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شکل ۱- موقعیت حوضه آبریز مسیل شاهین‌شهر در استان اصفهان و کشور 


۲-۲ روش پژوهش 

به‌منظور انجام مطالعات تعیین حدود بستر و تهیه مدل رقومی ارتفاع منطقه از نقشه ۰ (شرکت مهندسی 
مشاور آب پردازان نواندیش) که محدوده رودخانه و بخشی از دشت سیلابی را پوشش می‌دهد استفاده شده است 
I)‏ ۲)؛ و مدل توپوگرافی از رودخانه در قالب شبکه نامنظم مثلئی (TIN)‏ تهیه شد. تهیه داده‌ها در نرم‌افزار GIS‏ و 
انتقال آن‌ها به نرم‌افزار HEC-RAS‏ (ورژن 6.0.3 انجام گردید در مدل ریاضی HEC-RAS‏ روند اصلی محاسبات بر 
al‏ حل معادله انرژی یک بعدی استوار می‌باشد. افت‌های انرژی بر اساس افت اصطکاکی معادله مانینگ واگرایی و 


همگرایی " ضریبی که در تغییرات بار سرعت اعمال می‌شود ارزیابی می‌شود. در محبط GIS‏ با استفاده از ابزار HEC-‏ 


1 expansion 
2 Contraction 


AY‏ جغرافیا و مخاطرات محیطی شمارة اول 


Geo RAS‏ از مدل توپوگرافی رودخانه اطلاعات مناسبی چون مقاطع عرضی» موقعیت سواحل و فواصل بین مقاطع را 
با فرمت قابل‌استفاده در نرم‌افزار HEC - RAS‏ استخراج گردید. در اولین مرحله در HEC-Geo RAS‏ موقعیت 
سواحل رودخانه در مقاطع عرضی نیز در پیش‌پردازش داده‌ها در قالب لایه سواحل رودخانه تهیه و به مدل تعریف 
گردید در ادامه aU‏ ابعاد مسیر جریان به‌موازات لایه سواحل تهیه گردید. HEC-Geo RAS‏ تمام اطلاعات ژئومتری 
رودخانه با کمک L) RAS aU‏ فعال شدن (editor‏ تهیه می‌شود. لایه‌های RAS‏ ویژگی‌های ژئومتری زمینی را با 
رقومی کردن لایه‌های مختلف و سپس بازیابی دوباره دیتاها جمع‌بندی می‌کند Cual s)‏ ۲۰۱۲). در ادامه اطلاعات به 
HEC-RAS Low‏ منتقل شدند. یکی از پارامترها در مدل‌سازی» ضریب مانینگ است مقدار جریان در آبراهه‌ها تابعی 
از زبری جدار و بستر رودخانه می‌باشد این ضریب نمایانگر مقاومتی است که زبری جداره‌های رودخانه‌ها و مصالح 
آن در مقابل جریان و تغییرات آن ایجاد می‌نماید مقدار ضریب زبری با تأثیرپذیری از شرایط هیدرولیک دانه‌بندی مواد 
بستر» شکل مسیر و مورفولوژی رودخانه و همچنین عوارض و عوامل طبیعی و انسانی موجود در بستر رودخانه و 
سیلاب دشت‌ها در قسمت‌های مختلف آپراهه ممکن است متفاوت باشد (کوان ۰ 01( .جهت محاسبه ضریب زبری 
مانینگ از روش کوان(۱۹۵7) استفاده شده است ( جدول ۱). در مرحله بعد جهت تعیین دبی در دوره بازگشست‌های 
مختلف مقدار آن برای هر یک از زیر حوضه‌های شاهین‌شهر با استفاده از نرم‌افزار HEC-HMS‏ محاسبه شد (جدول 
(E‏ همچنین با توجه به شرایط موجود مبنای هیدرولیکی در محاسبه بستر این رودخانه جریان ماندگار در نظر گرفته 
شده است بر مبنای محاسبات انجام گرفته رژیم جریان رودخانه در شبیه‌سازی سیلاب اغلب مخلوطی زیر بحرانی و 
در موارد نادر فوق بحرانی می‌باشد؛ که با توجه به شرایط رودخانه مبنای محاسبه رژیم جریان رودخانه شاهین‌شهر 
رژیم مختلط در نظر گرفته شده است؛ و نتایج حاصل HEC-RAS‏ به محیط Civil3D‏ و سپس به محیط GIS‏ منتقل 


شده است و پهنه سیلاب در طول مسیل شاهین‌شهر بر اساس دوره باز گشت مختلف محاسبه گردید. 


جدول ۱- مقدار ضریب زبری مانینگ در مقاطع مسیل شاهین‌شهر 


ساحل راست بر ساحل چپ شماره مقاطع 
٩۰۷۹ - ۳ 4 ۷ HY‏ 
ío - 00 ۰/6۵ NY ۰/۳۵‏ 4" 
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شکل ۲- بخشی از TIN‏ مسیل شاهین‌شهر و مقاطع عرضی در محیط GIS‏ 


۱-۲-۲- محاسبه دبی در نرم‌افزار HEC-HMS‏ 

جهت آماده کردن اطلاعات ورودی به مدل 11۳0-115 ابتدا با استفاده از نقشه توپوگرافی حوضه در مقیاس 
۰ سازمان نقشه‌برداری تهیه شد و سپس بر اساس موقعیت حوضه نسبت به شناسایی ایستگاه‌های باران‌سنجی 
داحل و اطراف حوضه اقدام شد. با استفاده از نقشه‌های رقومی به مقیاس ۱:۲۵۰۰۰ مرز هیدرولوژیکی حوضه آبریز 
مسیل شاهین‌شهر تعیین و شبکه آبراهه‌های آن استخراج گردید که درمجموع حوضه موردمطالعه به ٤‏ زیر حوضه 
تقسیم‌بندی شد تمامی محاسبات در نرم‌افزار GIS‏ محاسبه گردید. شکل A)‏ شماتیک محیط نرم‌افزار HEC- HMS‏ 
جهت برآورد حداکثر دبی سیلابی با دوره برگست مختلف برای هر یک از زیرحوضهها را نان می‌دهد. در این 
نرم‌افزار جهت تحلیل و ایجاد هیدروگراف واحد برای هر یک از دوره بازگشست‌های مذکور ابتدا با وارد کردن حوضه 
آبریز منطقه موردمطالعه (مسیل شاهین‌شهر) به داخل محیط نرم‌افزار HEC- HMS‏ و سپس مشخصات فیزیوگرافی 
حوضه آبریز وارد گردید. بعد از تحلیلء نتایج خروجی حداکثر دبی سیلابی با دوره بازگشت‌های مختلف برای هر 
یک از زیرحوضه‌های حوضه آبریز موردهطالعه محاسبه شده است. 


HEC-HMS اطلاعات مدل در نرم‌افزار‎ ags 


به‌طورکلی در bE agi‏ ورودی دو سری داده برای مدل می‌بایست تهیه شود که عبارتنداز: 


Lal > At‏ و مخاطرات محیطی شمارة اول 


الف) مدل حوضه" در این روش جهت برآورد میزان دبی» از مدل‌سازی 552 با روش SCS curve number‏ 
استفاده گردید و مدل‌سازی بارش به رواناب نیز با روش (SCS unit hydrograph)‏ انجام شده است برای محاسبه 


pre‏ تفوذ در روش SCS-CN‏ ميزان cub‏ با استفاده از روابط زیر تخمین زده می‌شود: 


رابطه (۱) la = 02 * S‏ 
(S)‏ نگهداشت آب در سطح حوضه 
25400 
رابطه s == - 5 (Y)‏ 


شماره منحنی خاک (CN)‏ 

این پارامتر که بیانگر حصوصیات نفوذپذیری خاک حوضه می‌باشد polis‏ شماره منحنی خاک (CN)‏ با توجه به 
نوع کاربری اراضی و ویژگی‌های SE‏ تعیین می‌شود. برای تهیه فاکتور شسماره منحنی خاک نیاز به داده‌های نقشه 
گروه‌های هیدرولوژیک خاک کاربری اراضی. شیب منطقه در سطح حوضه آبریز می‌باشد. نقشه گروه‌های 
هیدرولوژیک خاک از داده‌های اداره کل منابع طبیعی استان اصفهان استفاده شده است؛ و نقشه کاربری اراضی با 
استفاده از تصویر ستتینل ۲ در محیط ag GIS‏ شده است. با توجه به نقشه گروه هیدرولوژیک خاک حوضه آبریز 
شاهین‌شهر che‏ چهار گروه خاک ABCD)‏ می‌باشد و انواع کاربری‌ها شامل تیپ‌های مرتعی, تیپ‌های جنگلی؛ 
زراعت آبی و دیم باغات در قلمرو مطالعاتی ملاحظه شده است. در این مطالعه مقادیر LCN‏ تلفیق نقشه‌های کاربری 
اراضی و گروه‌های هیدرولوژیکی خاک در محیط GIS‏ محاسبه گردید و سپس با تعیین ضریب رواناب حوضه آبریز 
بر اساس نوع پوشش سطحی حوضه و شیب اراضی در هر بخش تعیین شده است در میزان (CN)‏ ضریب شد و 
ely ale (CN)‏ تمامی زیر ia po‏ محاسبه گردید شای duty sell clays ACI glad‏ گروه‌هاق 
هیدرولوژیک و مقدار CN‏ حوضه مطالعاتی را نشان می‌دهد. 

همان‌طور که در جدول )1( ملاحظه می‌شود زمان تمرکز با استفاده از روش‌های کربای» کرپیچ و سازمان 
هوانوردی امریکاء برانسبای و csl ple) SCS‏ ۱۳۸۸) محاسبه شد و روش سازمان حفاظت خاک (SCS ( IS pol‏ و 
روش برانسبای ویلیامز مقادیر مناسبی را نشان می‌دهند. دلیل این امر به شرایط کاربرد این روش‌ها بازمی‌گردد این 
روش‌ها bare‏ برای حوضه‌های بسیار کوچک کاربرد dila‏ این در Sle‏ است که روش a SES‏ برانسبای- ویلیامز 
برای حوضه‌های موردمطالعه اعداد مناسب‌تری را ارائه می‌نماید؛ بنابراین زمان تمرکز حوضه و زیرحوضه‌های فوق با 
استناد به روش برانسبای- ویلیامز (علیزاده, ۱۳۸۸) در نظر گرفته شده است. 


1 basin model 
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مدل هواشناسی حوضه gal‏ 

ازجمله پارامترهای مهم در زمینه مطالعات آب و هواشناسی یک منطقه محاسبه شدیدترین بارش‌های چند دقیقه 
تا چند سععته در دوره بازگشت‌های مختلف است. این مهم در تمامی برنامه‌ریزی‌های محیطی منطقه به‌ویژه در 
مطالعات مربوط به سیلاب در رودخانه‌ها بسیار حائز اهمیت است. برای محاسبه شدت بارش منطقه موردمطالعه 
پارامترهای مورد نیاز ميانگین سالانه بارندگی و میانگین حداکثر بارش‌های ۲۶ ساعته ایستگاه جعف رآباد برآورد گردید. 
میانگین Le‏ بارندگی که با استفاده از آمارهای موجود محاسبه شد و حداکثرهای بارش ۲۶ ساعته نیز با استفاده از 
بارش روزانه ایستگاه فوق‌الذکر طی سال‌های آماری (۱۳۹۸-۱۳۷۰) استخراج گردید. 

برای محاسبه شدت بارندگی در مقیاس زمان (ساعت و دقیقه) و برای دوره بازگشت‌های مختلف نمودارهای 
شدت -مدت بارندگی تهیه شد (شکل € برای این کار ابتدا لازم است» PIG‏ در منطقه محاسبه شده که برای این 
منظور از روابط همبستگی استفاده می‌شود. Uke‏ با استفاده از داده‌های باران‌نگار ایستگاه‌های ایران مطالعه شده و مقداں 
0 به‌صورت زیر محاسبه شده است: 

(Y (رابطه‎ 

pro = 2.26(Pog) (b uL) °72 

که در Pagh Ol‏ متوسط حداکثر بارش‌های ۲۶ ساعته و Pyear‏ متوسط بارندگی سالانه منطقه برحسب میلیمتر 
است. برای نقاط مختلف Oly!‏ در مدت‌های بین ۱۵ دقیقه تا ۲۰ ساعت فرمول زیر ارائه شده است: 

(É (رابطه‎ 


P = [0.4524 + 0.2471 In(T — 0.6)](0.3710 + 0.6184¢°*48*) p£? 
بر حسب میلیمتر می‌باشد (علیزاده‎ PEP برحسب ساعت و دوره باز گشست برحسب سال و‎ t که در رابطه فوق‎ 


Ob sl آمارها در دوره بازگشت‌های مختلف بر اساس معادله‌های ارائه شده برای شرایط آب و هوایی‎ adal در‎ YM 
گردید و با توجه به ویژگی‌های فیزیوگرافی حوضه دبی در‎ HEC-HMS محاسبه گردیده است. نتایج وارد نرم‌افزار‎ 
دوره بازگشت‌های مختلف برای تمامی زیر حوضه‌ها محاسبه گردید.‎ 


حدول ۲- مقادیر حداکثر بارندگی ۲۶ ساعته ایستگاه جعفرآباد با دوره آماری (۱۳۹۸-۱۳۷۰) برای دوره 


jb‏ گشت‌های مختلف 


جعفر آباد ۲سال ۵ سال ۰ سال ۵ سال goo‏ | ۱۰۰سال | ۲۰۰ سال ۰ سال ۰ سال 
۵ دقیقه VAY "a/£v "WAR ۳ N41 yarw | TTY | YMWM‏ ۸0/۱ 
۰ دقیقه ۱۳/۳۸ ۱۸/۹۱ ۳۳/۵ £A/AA £o/AY 784 ywy ۷۲ YAN‏ 
۰ دقیقه ۷/۳ ۱۱/۳۵ ۱۳/۹۹۵ YV/oV Y£/RY ۲۲/۰۶ | MM YY‏ ۳۹/۹۵ 
۰ دقیقه 0/0۷ ۸/۱ 1/41 ۱۳/۹۱ ۳۲ ۱2/۵۲ ۱۸/۳۱ YY/£o yew‏ 
۲ ساعت 01/£ 4۷ VAYA yar Y£/44 ۱۳/۵۳ ۹۱ ۱/۷ MOV‏ 


^ جغرافیا و مخاطرات محیطی شمارة اول 
جعفرآباد ۲ سال ۵ سال ۰ سال سال | ۵۰سال | ۱۰۰سال | ۲۰۰ سال ۰ سال ۰ سال 
٤‏ ساعت MEZAN ۹/۰ NEA ۷/۹ vy ۱۳۷ £N 1/۸٦‏ ۱۱/۲ 
1 ساعت ۳/۹ N/^E DAE? v/Y\ Vos 0/۸۰ 0/*^ £/W Y/Yo‏ 
V‏ ساعت 1/4۷ ۳/۲ MAY vi VM WY o/Yí £NY Y/AV‏ 
A‏ ساعت VIVA VYY VEN 0/0۷ £NY ARÎ yyy YAY \NA‏ 
۲ ساعت A ۶ 1۹ YY ۳/۳۷ ۳/۳6 MAR! VEN‏ ۵/۵ 
VA‏ ساعت ۱/۰ £/oY ۷ Ava) ۳/۳۹ ۳/۹۰ T/0 ۳/۳۹ AS‏ 
۶ ساعت YA Y/£Y ۳/۳ YNY vty YAY yNY ۱/۱ “AY‏ 
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شکل ۳- نمودار شدت — مدت بارندگی ایستگاه جعفرآباد برای دوره بازگشت‌های مختلف 


حدول ۳- مساحت و lac‏ حوضه اصلی و زیرحوضه‌های حوضه آبریز مسیل شاهین‌شهر 


شماره زیر حوضه زیر حوضه ZI)‏ | زیر حوضه Z2)‏ | زیر حوضه Z3)‏ | زیر حوضه (ZA)‏ | کل حوضه 
مساحت oW£E/A YANA MÁÉYN ۸۷۱ YAYA (km2)‏ 
شیب ناخالص % UM VAY ZA vw WY‏ 
VYN VYN ۷۳/۵ wir ۷۳/۵ CN‏ 
طول آبراهه اصلی YAN 00/۸ WY ov" (km)‏ ۱۳۱/۵ 
کربای 1 vw VA‏ 0/۲ ۱/4 
YA £N vé VY 0/4 ES‏ 
o om pow ow pex] ese‏ ان .5 
w 1 ۷/۳ V/A "VA FAA‏ 
Y£/o WA ۱3/6۵ ۳/6۵ ۱۳/۳ SCS‏ 
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شکل ۶- نقشه کاربری اراضی حوضه آبریز مسیل شاهین‌شهر شکله- نقشه شیب حوضه آبریز مسیل شاهین‌شهر 
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شکل - گروه‌های هیدرولوژیک خاک حوضه شکل V‏ - تعیین مقدار (CN)‏ با احتساب 


آبریز مسیل شاهین‌شهر میانگین گیری وزنی 


ews 


ww ia M lee 
v= 


شکل ۸- شماتیک تعیین حداکثر دبی عبوری مسیل شاهین‌شهر بر اساس حداکثر دبی خروجی هر یک از 


زیرحوضه‌ها 


Lal > M‏ و مخاطرات محیطی شمارة اول 


جدول ۶- نتایج خروجی حداکثر دبی سیلابی (مترمکعب بر ثانیه) به ازای دوره بازگشت‌های مختلف برای کل 


حوضه ppl‏ شاهین‌شهر با استفاده از روش HEC-HMS (SCS)‏ 


دوره بازگشت (سال) 
روش‌ها 
EE o Yo E ۵ ۲‏ 
AVA- YAY yaw ۱۱/۳ HEC-HMS (SCS)‏ ۱۸/۵ ۱۳۳ 
-Y‏ نتایج و بحت 


با استفاده از سیستم اطلاعات جغرافیایی (TIN)‏ منطقه ساخته شده است و تمامی پارامترهای هیدرولوژیکی 
(مساحت. طول آبراهه اصلی CN‏ زمان تمرکز زمان تأخیر) در محیط GIS‏ محاسبه شد و تمامی این محاسبات وارد 
نرم‌افزار HEC-HMS‏ گردید دبی اوج لحظه‌ای برای ٤‏ زیر حوضه مسیل شاهین‌شهر محاسبه شد نتایج حاصل از 
مدل هیدرولوژی 11۳0-11۷15 در جدول )£( osla OLE‏ شده است و دبی اوج برای دوره بازگشت Yo‏ ساله AVA”‏ 
برای دوره بازگشت ۵۰ سال ۱۰۸/۵ و برای دوره بازگشت ۱۰۰ سال ۱۲۳/۱ مترمکعب بر A‏ محاسبه شده است؛ و 
تمامی این اطلاعات به‌عنوان ورودی نرم‌افزار HEC-RAS‏ برای مدل‌سازی عمق و سرعت آب و پهنه سیلاب در طول 
مسیل شاهین‌شهر استفاده شد. در محیط HEC-RAS‏ طول مسیر رودخانه شاهین‌شهر با ۲۵۰ مقطع که از DEM‏ 
استخراج گردید شبیه‌سازی شده است فاصله مقاطع به‌طور متوسط ۳۰ متر و عرض سطح مقاطع ۲۵۰ متر است و 
طول کل مسیر ۸ کیلومتر بوده است؛ و میزان ضریب زبری مانینگ برای کل مسیر رودخانه با توجه به ویژگی‌های 
سواحل سمت چپ و راست که اغلب دارای کاربری زمین بایره کشاورزی و مسکونی می‌باشد از حدود ۰/۱۷ تا 
۳ برای بستر و ۰/۲۶ تا ۰/۲۵ برای سیلاب‌دشت محاسبه شده است؛ و شرایط مرزی برای کل مسیر به‌صورت 
نرمال انتخاب گردید و نتایج حاصل از شبیه‌سازی برای سه دوره بازگشت YO‏ ۵۰ و ۱۰۰ ساله به‌صورت جریان 
حالت پایدار محاسبه شد و مسیل شاهین‌شهر با توجه به ویژگی‌های ژئومورفولوژیکی به ۳ بازه تقسیم‌بندی گردید. 

در بازه اول به طول ۳۲۰۰ متر بوده است میزان ضریب خمید گی این بخش امحاسبه شده است و با توجه به 
تقسیم‌بندی بریلی و فریس (۲۰۰۵) جزء الگوی مستقیم می‌باشد (حسین زاده و همکاران (YE‏ جنس بستر 
رودخانه ماسه و سیلت تشکیل شده است در ابتدای بازه دارای کاربری بایر و سپس کاربری مسکونی می‌شود و 
متوسط شیب بستر ۰/۲ درصد بوده است پهنه سیلاب در بازه اول در دوره بازگشت ۵ JU.‏ از حد مقاطع فراتر نرفته 
است اما در دوره بازگشت ۵۰ ساله ۱۰ تا ۲۰ متر و در دوره بازگشت ۱۰۰ سال متوسط ۳۰ متر از حد مقاطع فراتر 
رفته است در بازه اول با توجه به نقشه کاربری تهیه شده حاشیه مسیل شاهین‌شهر غالب منطقه دارای کاربری 
مسکونی ory‏ است که با توجه به پهنه گسترش سیلاب لزوم برنامه‌ریزی بیشتری را در این زمینه می‌طلبد. شکل )4( 
(اشکال call‏ و ب) 


سال یازدهم تحلیل دینامیکی مخاطره سیلاب با استفاده از ... ^ 


در بازه دوم طول مسیل شاهین‌شهر ۳ کیلومتر بوده است و میزان ضریب خمیدگی OF‏ ۱/۰۹ محاسبه شده است 
که ansa OSI‏ کی یادا کین افو gs: A aa‏ حاقیه سیل کارا کاریری cdi‏ بره 
است متوسط شیب بستر در این بازه ۰/۱ محاسبه شده است متوسط عمق مقطع در این بازه ۶ متر بوده است پهنه 
سیلاب در دوره بازگشت‌های مختلف در این بازه گسترش کمتری داشته است و در دوره بازگشت ۲۵ ساله پهنه 
گسترش سیلاب تا حدود ۱۹۰ متر از حد مقاطع فراتر رفته است. شکل(۱۰» (اشکال ۱۰ الف و ب) 


تصویر گوگل ارث عسیل شاهین شهر در بازه اول 
راهنمای نقشه 


بهنه سیلاب در دوره بازگشت ۲۵ ماه 5 
پهنه سیلاب در دوره باز گشت ۵۰ ساله 


پهنه سیلاب در دوره باز گشت ۱۰۰ ساله 


£ |025 05 O75 1 
;jj- -—— ——88 — 


5 ~ = = ona ate, و"‎ È or 
er = pei m: -[— ate 
— 3 ath BS cte = no ee). eS. Tn P - A 


شکل -٩‏ نقشه بهنه‌بندی سیلاب مسیل شاهین‌شهر در دوره بازگشت ۲۵ و ۵۰ و ۱۰۰ساله 


.4 جغرافیا و مخاطرات محبطی شمارة اول 


شکل ٩‏ الف و ب- مقطع عرضی مسیل شاهین‌شهر در دوره بازگشت 0« OF‏ و ۰ سال و پهنه سیلاب در دوره 
بازگشت cae YO‏ ١۵د‏ »بات 
تصوبر گوگل ارث مسیل شاهین شهر در بازه دوم 


پهنه سیلاب در دوره بازگشت ۵۰ ساله 


بهنه سیلاب در دوره باز گشت ۱۰۰ ساله 
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زگشت ۲۵ و ۵۰۰ و JU‏ 


شکل ۱۰- نقشه پهنه‌بندی سیلاب مسیل شاهین‌شهر در دوره با 


سال ea‏ تحلیل دینامیکی مخاطره سیلاب با استفاده از .... ٩۱‏ 


شکل ۰ الف و ب- مقطع عرضی مسیل شاهین‌شهر در دوره بازگشت 0« ۰ و ۰ سال و پهنه سیلاب در 


دوره بازگشت YO‏ مت conte‏ ۸ات 


در بازه سوم طول بازه ۲ کیلومتر بوده است مسیل شاهین‌شهر در این بازه از میان زمین‌های بایر و کشاورزی در 
m S‏ کب وواه ور این سک ۷/۱ IERI TET‏ رکه امش کیب 
و کاهش سرعت Whe‏ پهنه سیلاب به حداکثر گسترش خود نسبت به قسمت‌های قبلی رسیده است. در این 
قسمت پهنه سیلاب در دوره بازگشت ۲۵ ساله» به حداکثر گسترش خود حدود ۱۰ تا ۰ متر از حد مقاطع گسترش 
پیدا کرده است در دوره بازگشت ۵۰ سال از ۱۲ تا ۱۸۰ متر و در دوره بازگشت ۱۰۰ ساله تا ۲۰۰ متر از حد مقاطع 
سیلاب گسترده شده است: 

ا و اک نی روات [Lo‏ قافن ر مات و ala‏ ل et‏ ارق SUS‏ ر 
فرسایش بوده است خطر ناشی از سرعت‌بالا و انتقال رسوب در انتهای بازه موردمطالعه با توجه به کاهش شیب 
مخاطره‌ای که ممکن است ادیده گرفته شود انباشت رسوبات (با توجه به جنس بستر مسیل) می‌باشد که تحقیقات 
بیشتری Lb‏ در مورد حمل‌ونقل رسوبات با استفاده از مدل‌های عددی جهت کاهش مخاطرات انجام شود. 
شکل(۱۱» ca SIKH)‏ و ب) 


ay‏ جغرافیا و مخاطرات محیطی Sess‏ ول 


تصوبر گوگل ارث مسیل شاهین شهر در بازه سوم 
ر اهنمای نقشه 


بهنه سیلاب در دوره بازگشت ۲۵ اه BE‏ 
پهنه سیلاب در دوره باز گشت ۵۰ AIL,‏ ۳ 


پهنه سیلاب در دوره بازگشت ۱۰۰ ساله 
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شکل ۱ الف و ب- مقطع عرضی مسیل شاهین‌شهر در دوره بازگشت 0۵ ۰ و ۰ سال و پهنه سیلاب در 
دوره بازگشت ——\rs ,--0+ ,—— YO‏ 


سال یازدهم تحلیل دینامیکی مخاطره سیلاب با استفاده از .... ay‏ 


در کل بازه موردمطالعه متوسط عمق آب در دوره بازگشت ۲۵ ساله ۱/۸۸ در دوره بازگشت ۵۰ ساله VAY‏ در 
دوره بازگشت ۱۰۰ ساله ۲/۰۵ متر رسیده است و علاوه بر آن حداکثر عمق آب در دوره بازگشت ۱۰۰ سال به میزان 
۷ متر در بازه سوم رخ داده است. در کل بازه موردمطالعه مسیر رودخانه ساماندهی شده در انتهای بازه موردمطالعه 
عمق آب با توجه به دوره بازگشت‌ها و میزان دبی افزایش پیدا کرده است و پهنه گسترش سیلاب وسیم‌تر از سایر 
قسمت‌ها شده است؛ که در جداول (۵الف و ۵ب) متوسط شیب» عمق. مساحت مقطع عرضی بیان شده است. 


جدول ۵ الف - مقایسه میانگین مساحت مقطع عرضی و عرض مقاطع در بازه‌های مختلف مسیل شاهین‌شهر 


دوره بازگشت ۵ سال 
متوسط متوسط متوسط عمق " 
عرض دبی عمق اب G9)‏ 
نام بازه مساحت مقطع | شیب بستر مقاطع عرضی 
G2) uU‏ متوسط 
عرضی (متر) (درصد) (متر) متوسط حداکثر 
کل باز 
بازه اول ۱۱/۸۷ ۰/۲ ۱۸/۹ ۱/۳۳۸ YIM YAY‏ 
بازه دوم ۱۳/۹ ۰/۱ )14/4 Y/AY ۱/۸۱ (AE‏ ۱/۸۸ 
بازه سوم ۱۳/۸۹ ۰/۱ ۱/۸۲۳۷ ۰/4 ۱/۸ 01/£ 
جدول ۵ ب- مقایسه میانگین عمق آب در بازه‌های مختلف مسیل شاهین شهر 
دوره بازگشت ۰ سال دوره بازگشت ۱۰۰ سال 
عمق آب (متر) عمق آب (متر) 
نام بازه 
متوسط کل متوسط کل 
متوسط حداکثر متوسط حداکثر 
بازه بازه 
بازه اول £N Y/YA‏ ۳/۳ 1۹ 
بازه دوم Y/AA VAY‏ ۱/۹ ۲ ۳/۳۹ ۲/۰0 
بازه سوم ۱/۷۵ £N ۱/۷۹ LAE‏ 


“Wyte‏ مساحت پهنه سیلاب در دره بازگشت‌های مختلف مسیل شاهین‌شهر 


دوره بازگشت مساحت پهنه سیلاب به (کیلومتر مربع) 
Yo‏ ۰/۹۸ 
vir 0۰‏ 


۱/۳ Um 


Af‏ جغرافیا و مخاطرات محیطی شمارة اول 


ct‏ نتیجه گیری 

sib وان‎ pM ie alte Sealy ull ae co ee کی از‎ Lea bu 
ایران مستعد مخاطره سیلاب است و‎ c gl می‌کند. به دلیل شرایط توپوگرافی و آب‌وهوای و ضعف مدیریت حوضه‎ 
در سال‌های اخیر حوادث شدید سیل را تجربه کرده است اولین گام اساسی در مدیریت سیلاب تولید نقشه‌های‎ 
زیر حوضه‎ ٤ مناطق مستعد سیلاب بوده است. حوضه آبریز مسیل شاهین‌شهر با توجه به شرایط توپوگرافی به‎ 
تقسیم‌بندی شد و سپس زمان تمرکز برای هر یک از زیر حوضه‌ها محاسبه شد و زمان تمرکز به دست آمده به‎ 
روش‌های کربای. کرپیچ و سازمان هوانوردی امریکا بسیار کوچک‌تر از مقادیر حاصل از روش سازمان حفاظت خاک‎ 
امریکا و روش برانسبای ویلیامز می‌باشد و زمان تمرکز حوضه با توجه به وسعت حوضه و نتایج به ست آمده برای‎ 
۱۳/۲۱ هریک از زیرحوضه‌های فوق با استناد به روش برانسبای- ویلیامز در نظر گرفته شد که برای کل حوضه‎ 
حداکثر دبی‎ HEC-HMS(SCS) J ile 5 و با استفاده از‎ SCS ساعت محاسبه شده است و سپس با کمک روش‎ 
گردید. مسیل‎ HECRAS سیلابی برای دوره بازگشست‌های ۲۵ ۵۰ و ۱۰۰ ساله محاسبه شد و وارد نرم‌افزار‎ 
پهنه‌بندی‎ HEC-RAS شاهین‌شسهر با ضریب خمیدگی کل ۱/۰۹ دارای الگوی سینوسی می‌باشد. با اجرای مدل‎ 
ساله ۵۰ و ۱۰۰ ساله اجرا شد که متوسط پهنه گسترش سیلاب در دوره‎ Yo سیلاب برای سه دوره بازگشت‎ 
بازگشت‌های ۲۵ ۵۰ و ۱۰۰ ساله از ۲ تا ۲۰۰ متر از حد مقاطع گسترش پیدا کرده است در بازه سوم با توجه به‎ 
اینکه مساحت مقطع عرضی بیشتر از سایر بازه بوده است ولی با کاهش شیب پهنه گسترش سیلاب در این بازه از‎ 
سایر بازه‌های گسترده‌تر شده است که با توجه نوع کاربری در این بخش که اغلب کشاورزی و مسکونی بوده است‎ 
پهنه وسیعی تحت مخاطره سیلاب در دوره بازگشت‌های مختلف قرار می‌گیرد. اگر محدوده پهنه سیلاب در دوره‎ 
بازگشت‌های مختلف را جز منطقه آسیب‌پذیر بدانیم اغلب مناطق مسکونی در معرض مخاطره سیلاب قرار گرفته‌اند‎ 
و به برنامه‌ریزی بیشتری نیاز دارد این مطالعه با نتایج سایر محققان ازجمله عز ( ۲۰۱۸) و شسفیعی مطلق و عبادتی‎ 
plat do عطافت‎ (La Raw بتافی و‎ Ob) ail Rew d ae CU) a Rae یضار و‎ CR) 
پیش‌بینی شده است و همچنین‎ HEC-RAS مطالعات بیان می‌کنند که مناطق مستعد سیلاب به‌وسیله مدل عددی‎ 
در برآورد دبی در حوضه‌هایی که فاقد ایستگاه آبسنجی بوده است‎ HEC-HMS مطالعات نشان می‌دهد که مدل‎ 
مناسب می‌باشسد. علاوه بر این» شبیه‌سازی» اطلاعات بسیار مفیدی مانند ارتفا سسرعت و زمان جریان برای‎ 
تصمیم گیرندگان فراهم می‌کند. این نتایج می‌تواند به تولید نقشه‌های مناطق مستعد سیل منجر شود این اقدامات‎ 
سیلاب را از طریق پیش‌بینی خطی میزان آسیب‌پذیری ناشی از سیل را کاهش دهد. شناسایی‎ SE مدیریتی می‌تواند‎ 
مناطق مستعد سیل می‌تواند جوامع تهدید شده و تصمیم گیرندگان را در سطوح مختلف برای جلوگیری از مخاطره‎ 
آگاه کند و بنابراین این امر کلید مدیریت خطر سیل است.‎ 


سال یازدهم تحلیل دینامیکی مخاطره سیلاب با استفاده از .... qo‏ 
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